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Graphene quantum dots are considered as promising alternatives to quantum dots in III-V semiconductors, 
e.g., for the use as spin qubits due to their consistency made of light atoms including spin-free nuclei which both 
imply relatively long spin decoherene times. However, this potential has not been realized in experiments so far, 
most likely, due to a missing control of the edge conhgurations of the quantum dots. Thus, a more fundamental 
investigation of Graphene quantum dots appears to be necessary including a full control of the wave function 
properties most favorably during transport spectroscopy measurements. 

Here, we review the recent success in mapping wave functions of graphene quantum dots supported by metals, 
in particular Ir(lll), and show how the goal of probing such wave functions on insulating supports during 
transport spectroscopy might be achieved. 


Graphene has moved in short time from first preparation 
as a small flake Bl towards applications, partly already 
realized and partly anticipated for the near future, such as 
high frequency transistors O, supercapacitors |i6j], or touch 
screens m as recently reviewed by Ferrari et al. (H. Most of 
the applications are based on the exceptional properties of the 
material including true two-dimensionality, inertness, high 
room temperature mobility, large thermal conductivity, and 
extreme mechanical breaking strength ||9l . 

Another exciting plan is to use graphene quantum dots (QDs) 
as spin qubits ca. The basic requirement is a very long spin 
coherence time im, which might exist in graphene ca due 
to the absence of hyperfine coupling in isotopically pure ma¬ 
terial and the small spin-orbit coupling ca. First graphene 
QDs have been produced early on and probed by transport 
measurements ifTH - IT^ . meanwhile being able to determine 
charge relaxation times of about 60-100 ns mi. However, 
since graphene provides no natural gap, it is difficult to 
control the electron number lIT^ . Moreover, the 2D sublattice 
symmetry makes the QD properties very susceptible to the 
atomic edge configuration cni in contrast to the behavior 
in conventional QDs made, e.g., out of GaAs. As a result, 
chaotic Dirac billiards have been predicted for graphene 
QDs |[T^ and were even claimed to be realized IIHlJOl, i.e. 
the wave functions are probably rather disordered. To get 
more control on graphene QDs, the QD edges must be well 
defined and a more fundamental understanding of the QD 
properties is mandatory. More recent transport investigations 
on graphene nanoconstrictions indeed support the idea of 
a decisive role of the edges for the localization properties 
within graphene nanostructures ED. 

A direct insight into QD properties is provided by scanning 
tunneling spectroscopy (STS). STS can map the squared 
wave functions of QDs ca and, at the same time, determines 
the shape of the QD atom by atom. STS has also been applied 
successfully to graphene samples EH- However, conven¬ 
tional STS is restricted to conducting surfaces, such that first 


investigations of QD wave functions have been performed 
on graphene nanoislands prepared in Ultra-High-Vacuum 
(UHV) on metals El 12414291 . More peculiar quantum dots 
have also been probed by scanning tunneling microscopy 
(STM) in UHV such as the ones induced by tip-induced strain 
on suspended graphene areas EOl , the ones being present 
within localized areas of a quantum Hall sample ED, or the 
ones confined in the suspended areas of the strongly buckled 
graphene on Ru(OOOl) 02 • 

Here, we firstly review the results of wave function mapping 
of graphene QDs on Ir(lll), which includes the probing of 
the zig-zag edges, appearing to be the energetically favorable 
edge type for this substrate ES- It is found that the graphene 
QDs exhibit squared wave functions El 12414261 . which are a 
mixture of graphene states from the Dirac cone located in a 
projected band gap of the Ir(lll) E4l and hole-type Ir(lll) 
surface states located around F and exhibiting a strong 
Rashba-type spin splitting E5l . The strength of the two 
contributions within the QD can be disentangled by Fourier 
transformation of the real space data, if atomic resolution is 
achieved ED- The contribution from the graphene Dirac cone 
state can, moreover, be enhanced by oxygen intercalation 
between the graphene nano-islands and the Ir(lll) surface 

ED- 

Interestingly, the edge configuration of the graphene quantum 
dots on Ir(lll) is quite unique by the hybridization of the 
graphene Pz orbitals with Ir dz^^ orbitals, the latter being 
located rather exactly at the Fermi level This, on the 
one hand side, suppresses the appearance of the famous edge 
state at zig-zag edges E6l and, on the other hand side, leads 
to a soft confinement of the graphene quantum dot states El. 
Thus, the confinement is not atomically abrupt, as one might 
anticipate from the abrupt end of carbon atoms at the rim 
of the island, but proceeds over 4-5 lattice constants. As a 
result, the interference effects, coming from the coupling of 
the two valleys K and K’ and typically making the quantum 
dot wave functions chaotic, are suppressed EEQl, which. 
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in turn, makes the wave functions rather regular. Indeed, 
calculations of the atomically identical quantum dots, which 
are freely suspended and, thus, not prone to the interaction 
with the substrate, show strongly irregular wave functions 
distinct from the experiment. Only the inclusion of the soft 
confinement potential could recover the experimental results. 
Thus, it is experimentally demonstrated that soft confinement 
leads to a significant regularization of states in graphene 
quantum dots and, thus, should be part of a strategy to exploit 
the potential of graphene quantum dots, e.g., with respect to 
spin qubits nni. 

In addition, the soft confinement is able to suppress the 
graphene edge states. Such states with almost fiat bands at the 
Dirac point energy are an intrinsic property of graphene 
zigzag edges, at least, if unsupported and unreconstructed 
E3iQl- Experimental evidence for such edge states has 
indeed been found on different substrates such as HOPG 
EB ga, Si(lOO) El, Au(lll) El 03, Cu(OOl), if the 
graphene is laterally interfaced with BN, 1461 . arguably 
Pt(lll) E3, SiC(OOOl) EH, and Si02 ESlEOl. Spectro¬ 
scopic features at the step edges of Ir(lll) covered with 
graphene have also been interpreted in terms of graphene 
edge states ED, albeit it is known that the graphene layer is 
not interrupted at the Ir step edges EH. 

Density functional theory calculations indeed show that the 
edge state is present on Cu(lll), Ag(lll), and Au(lll) 
1531, but not on Ir(lll) |[36l|54l- It is, moreover, shown that 
H-terminated zig-zag edges on Au(lll) are ferromagnetic, 
i.e. the edge state is located at ES Indirect evidence 
for such ferromagnetism of the edge state has been published 
omi, but without chemical control of the edge termination 
so far, — being possible in principle 1551 . Such edge state 
magnetism might be exploited for a multitude of spintronic 
applications 156U591 , but, on the other hand, it is detrimental 
for quantum dot functionality. Thus, the edge state suppres¬ 
sion by soft confinement, as on Ir(lll), again improves the 
reliability of the quantum dot wave functions. 

Interestingly, the edge state can be recovered by laterally 
interfacing the graphene zig-zag edges on Ir(lll) with BN 
and subsequently intercalating Au between the graphene/BN 
layer and the Ir(l 11) surface l6Ql . 

Albeit these STS studies of graphene QDs on metals open 
the door towards a better understanding of QD wave function 
properties, they are still relatively remote from the transport 
investigations dElED. A direct imaging of the wave 
functions responsible for the transport signatures would be 
more favorable. This goal is also pursued by scanning gate 
microscopy ED, but the achieved spatial resolution so far 
only allows to locate the center of the corresponding wave 
functions without any details EDisa. Combining transport 
and scanning tunneling microscopy is possible in principle, 
if the nanostructures are produced without resist as, e.g., by 
anodic oxidation using an atomic force microscope (AFM) 
||63]| . This technique has been shown to work on graphene, 
too 1641 . and to produce operating quantum dot structures 
lf65l . Besides cuts, where the graphene is removed l(64l . the 
graphene can also be oxidized 1661 or hydrogenated 1671 . 
which eventually could lead to smaller, insulating device 
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FIG. 1. (a) Sketch of exposure of C 2 H 4 to Ir(lll) at 300 K; (b) 
sketch of subsequent heating in order to produce graphene; (c) 3D 
representation of a (100 x 100) nm^ STM image of Ir(lll) covered 
by monolayer graphene islands; V = —0.3 V, / = 0.3 nA; upper 
inset: atomically resolved STM image of graphene island with zoom 
into zigzag edge having a sketch of the graphene lattice overlaid, 
V — 0.7 V, / = 20 nA; lower inset: zoom into Ir(lll) with atomic 
resolution, E = 0.9 V, / = 0.5 nA; (d) ARPES spectrum of clean 
Ir(lll) recorded along TM direction with positions of the K points 
of iridium and graphene marked as Kir and Kg, respectively; Si — 
S3 are surface states; (e) ARPES spectrum of Ir(lll) covered by 
graphene along the same direction as in (a); Dirac velocity ud, Dirac 
cone replica (R) and minigaps at the crossing points of original and 
replicated Dirac cones (arrows) are marked. |[3l[34l ((c), (d) courtesy 
of C. Busse, University of Cologne) 


boundaries. Importantly, this method does not leave any resist 
on the surface, which would disturb the STS measurements 
significantly due to a remote charging of the resist by the tip, 
which electrostatically would back-act on the area probed by 
STS (681. Combining this AFM based structuring method 
with a contacting method, which does not use any resist either 
as, e.g., microsoldering by In |[68]|69l, provides an adequate 
nanostructure to be probed both, by STS and transport 
spectroscopy simultaneously. 

In order to scan across the insulting barriers, the scanning 
tip has to be mounted on a tuning fork (TOl , such that the 
feedback can operate in scanning force mode while tunneling, 
or, alternatively, a capacitive feedback must be used ED 
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FIG. 2. (a)-(c) Calculated squared wave functions of the island probed in (e)-(h) using third-nearest neighbor tight-binding calculations which 
ignore the substrate; energies are marked; (d) calculated LDOS of the same island (consisting of 6 wave functions) at the energy marked; 
(e)-(h) dlIdV maps of the island at different energies E = V ■ e sls marked; / = 0.3 nA, Umod = 10 mV; question marks mark the bright 
rim observed in experiment, but not in the TB calculations; (i)-(j) calculated squared wave functions of the same island, which are not found 
in the experiment as marked by question marks; all calculations are with abrupt edge. (S) 


I. GRAPHENE QUANTUM DOTS ON IR(lll) 

A. Wave function mapping 

Monolayer graphene islands acting as QDs can be prepared 
by exposing clean Ir(lll) for 4 min to a pressure of 10“^ 
Pa of C 2 H 4 at 300 K and subsequent annealing to 1320 K 
(30 s) 117211 as sketched in Fig. [^a) and (b). The resulting 
islands (graphene QDs) have diameters of 2 — 40 nm and are 
completely enclosed by zig-zag edges as shown in Fig. [^c). 
Fortunately, Ir(lll) exhibits a projected band gap in the 
region of the graphene Dirac cone such that graphene 
QD states can be observed EiaEa. The band gap located 
between the surface states Si and S 2 is visible in the angular 
resolved photoelectron spectroscopy (ARPES) data shown 
in Fig. Bd)- After preparing a complete monolayer of 
graphene, the Dirac cone of graphene becomes visible within 
that band gap showing the well-known linear dispersion with 
a Fermi velocity of = (1.0 ± 0.05) • 10^ m/s 12911341 
(Fig. [2^))- The moire superstructure between graphene 
and Ir(lll) induces replica of the Dirac cone (marked by R) 
and corresponding minigaps at their crossing points (arrows) 
located at higher binding energy M- More importantly, 
the extrapolated Dirac point is at 0.1 eV above Ep l34l . 
overlapping with the Ir surface state Si, which has 
character |36l, thus, the Si orbitals point directly to the pz 
orbitals of graphene. Consequently, hybridization of Si with 
the Dirac cone and, thus, a gapping of the Dirac cone at Ed is 
likely and has indeed been found in density functional theory 
(DFT) calculations |3l[3l. Moreover, the projected band gap 
does only exist for the hole part of the band structure, while 
the electron part overlaps with bulk bands from Ir(l 11 ) lO. 


The local density of states (FDOS) of the QDs can be 
mapped by STS in constant current mode EElllBl. We 
use a STM operating at T = 6 K f73l and a lock-in tech- 
nique with modulation frequency u = 1.4 kHz and amplitude 
Vmod = 10 mV resulting in an energy resolution 6Ests ~ 
■ ksTf -t- (1.8 • eV^odY = 18meV El. This en- 
ergy resolution is much better than the natural peak width of 
the confined states of = 0.1 — 0.4 eV, which increases lin¬ 
early with \E — Ey\. For dl/dV curves, we stabilize the tip 
at sample voltage 14 tab and current /stab before switching off 
the feedback. Fig. [^e)-(h) show the FDOS patterns observed 
for a particular graphene QD in comparison with third-nearest 
neighbor tight binding (TB) calculations, which neglect any 
infiuence of the substrate (Fig. [^a)-(d)). The exact atomic 
configuration of the graphene QD deduced from the STM data 
has been used for the calculations. While single fourfold de¬ 
generate wave functions contribute at lower binding energy 
(Fig. I^a)-(c)), the FDOS at higher binding energy consists 
of six overlapping fourfold degenerate wave functions. A rea¬ 
sonable correspondence of the wave function symmetries be¬ 
tween TB results and STS results is found, but there are two 
decisive differences. 

Firstly, the experimental FDOS exhibits a bright band at the 
rim of the island (question marks in Fig. |^f)-(h)) which is not 
present in the TB data. Secondly, the strongly disordered wave 
functions expected from the sublattice symmetry breaking at 
the rim of the QD |[T8l[38l is found in the TB calculations (Fig. 
I^i) and (j)), but not in the STS data. 

The first discrepancy is related to the penetration of the sp- 
type Ir surface state into the graphene QD, while the sec¬ 
ond one is due to the soft confinement resulting from the hy¬ 
bridization of the graphene pz orbitals with the Ir(lll) dz 2 
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surface state. 



FIG. 3. (a) Energy gap /\E at the Dirac point versus graphene- 
Ir distance d as deduced from DFT calculations; upper right inset: 
band structure around for d = 0.3 nm with /S.E indicated, grey 
area corresponding to projected bulk bands of Ir and thick black lines 
mark graphene states; lower left inset: sketch of a graphene island 
(red) on Ir(l 11) (blue grey) with distances marked according to DFT 
calculations and X-ray standing wave data (331; (b)-(e) LDOS (= 
1^1^) deduced from TB calculations for individual confined states; 
(b), (c) E-Ey = -570 meV; (d), (e) E-Ey = -620 meV; (b), (d) 
calculated with abrupt edges; (c), (e) calculated with additional soft 
edge potential; (f)-(h) dl / dV images recorded at energies E — V -e 
as marked, 7 = 0.2 nA, ymod = 10 mV; (i)-(k) LDOS maps from TB 
calculations with soft edge potential at energies indicated, number 
of states contributing to the LDOS are marked. O ((a) courtesy 
of R. Mazzarello, RWTH Aachen University, (b)—(e) and (e)—(k) 
courtesy of F. Libisch, TU Vienna) 


with (Jz being the Pauli Matrix acting on the sublattice de¬ 
gree of freedom. Figure [^b)-(e) demonstrate that this on¬ 
site potential leading to soft confinement across about 5 lat¬ 
tice constants indeed removes the irregularities of the squared 
QD wave functions. Reasonable agreement between STS data 
and the TB data with soft confinement is achieved for smaller 
islands (Fig. |^f)-(k)), where the outer rim of Fig. [^f)-(h) ap¬ 
pears to be less strong. Notice that the moire potential, known 
from the size of the minigaps (see Fig. [Je)), is included as an 
on-site potential, too O. 

In larger islands as in Fig. Qa) and (b), the bright band at the 
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In order to take the hybridization into account, we per¬ 
formed density functional theory calculations of extended 
graphene sheets with different distance between the graphene 
layer and Ir(III). The resulting band structure (see inset of 
Fig. [^a)) exhibits a band gap /S.E at the Dirac point of the 
graphene related states. This /S.E increases with decreasing 
graphene-iridium distance d as displayed in Fig. [^a). The 
change of this distance at the rim of a graphene island has 
been deduced previously by favorable comparison of DFT cal¬ 
culations and x-ray standing wave experiments 13^ . It is lo¬ 
cated within the outer I nm of the graphene island and changes 
from 0.16 nm at the last atom to 0.34 nm in the interior of the 
graphene island as shown in the lower right inset of Fig. |^a). 
We take this into account by an on-site potential Vi at each 
lattice site i in the TB calculation reading 09): 


FIG. 4. (a)-(b) dl/dV maps of a graphene QD recorded at the volt¬ 
ages marked; 27 x 30 nm^, I = 0.5nA, Unod = 10 mV; blue line 
marks the orientation of the profile line shown in (c); (c) profile line 
along the line marked in (a) with deduced wave lengths Aout (Ain) 
outside (inside) the QD indicated; (d) resulting dispersion relations 
E{Ak = tt/A in/out) inside (red) and outside (blue) of the QDs as 
deduced from several QDs; full lines are linear fits with resulting vy> 
and relative energy offset marked; (e) ARPES data of Ir(lll) show¬ 
ing a Rashba split hole-type surface state, energy maximum (binding 
energy: 340 meV) is marked; the blue fit line from (d) (dashed line) is 
overlaid; (f) ARPES data of Ir(lll) covered with a complete mono- 
layer of graphene; energy maximum (binding energy: 190 meV) is 
marked; the red fit line from (d) (dashed line) is overlaid (35). O 
((e), (f) courtesy of A. Varykhalov, Helmholtz-Zentrum, Berlin) 

rim of the QDs develops into a standing wave at the rim, which 
is also apparent with similar wave length at the outer side of 
the island. The wave length dependence on energy E = e -V 
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is determined everywhere, where at least two maxima of the 
standing wave can be discriminated as shown in Fig. Be). 
The wave length of a standing wave A is related to the wave 
number k of the electronic state by k = tt/A giving access 
to the dispersion relation E{k) of the corresponding states in¬ 
side and outside of the graphene QD. It is shown in Fig. [^d) 
revealing the same slope inside and outside the island. As¬ 
suming a linear dispersion E = H • • k, one reveals a Dirac 

velocity = (4.5 ± 0.5) • 10^ m/s, which is less than half 
of the value known for the graphene layer on Ir(l 11) (see Fig. 
[Je)). Moreover, the crossing point of the fit lines in Fig. Qd) 
are at 0.15 — 0.3 eV below E^, while the extrapolated Dirac 
point of graphene is at 0.1 eV above Ep. Thus, the standing 
waves are not related to the Dirac cone of graphene. In con¬ 
trast, they fit nicely to an Ir(l 11) surface state of sp-type called 
So and being located around F, which is shown in Fig. Be) 
and (f) Ea. So shows a strong Rashba-type spin splitting and 
a nearly linear dispersion away from F. It survives the cov¬ 
erage of a graphene monolayer and is rigidly shifted by about 
0.15 eV upwards. It is well known that the standing waves 
observed in STS are not susceptible to the Rashba-type spin 
splitting l75l . such that the wave length of a standing wave is 
related to the average k of both bands. The averaged slope of 
both bands corresponds to = 5.5 • 10^ m/s and crossing 
points with F at 0.1 — 0.25 eV below E^ in good agreement 
with the standing wave data obtained by STS at the graphene 
QD rims. The relative energy shift of 0.15 eV upwards af¬ 
ter graphene coverage is also in agreement with the STS data. 
The fit results from the STS data are overlaid as dashed lines 
in Fig. I^e) and (f) substantiating the good agreement. Re¬ 
maining small differences are probably related to changes of 
the electrostatics at the rim of the islands with respect to the 
extended film probed by ARPES. 

In order to corroborate the conclusion that Sq contributes to 



FIG. 5. (a) Calculated LDOS of So for Ir(l 11) covered with graphene 
2 Lt E = — 0.4 eV along the direction perpendicular to the surface; 
the intensity in the Ir surface layer Iir and the intensity within the 
graphene layer Ic as used for the determination of = Ic/hr 
are marked; (b) relative intensity R = Ic/hr of So and S2 in 
graphene as deduced from STS data (squares) and from DFT cal¬ 
culations (So: circles, S 2 : triangles). O (Courtesy of R. Mazzarello, 
RWTH Aachen University) 

the LDOS patterns of the graphene QDs, we analyzed the ver¬ 
tical distribution of the So states as calculated by DFT. As 
shown in Fig. [^a), Sq indeed penetrates into the graphene 
layer with a few percent of its intensity. We compared the 


strength of Sq in the surface layer of Ir /jr and the strength of 
So in the graphene layer Iq giving a ratio Rbft = Ic/hr- 
This ratio displayed in Fig. [^b) is 10-15 %, while it is neg¬ 
ligible for other d-type surface states of Ir(lll) as, e.g., S 2 . 
The ratio has been compared with the relative amplitudes of 
the standing waves inside and outside the graphene QD (see 
Fig. Qc)) after recalibrating for different tip-surface distances 
above Ir(lll) and above graphene being deduced from the 
known height of graphene above Ir(lll) (0.34 nm. Fig. |^a)) 
and recorded I{z) curves. As shown in Fig. [^b), the agree¬ 
ment between STS and DFT data is very good, substantiating 
that So contributes to the LDOS patterns in the QD. 


B. Disentangling graphene Dirac cone states and Ir surface 
states 

The contributions from Sq and from the graphene Dirac 
cone to the LDOS of the graphene QDs depends on the 
details of the micro tip. Firstly, we show an example where 
the LDOS is dominated by the graphene states as for the 
images displayed in Fig. |^f)-(h). Figure |^a) shows the first 
two peak energies deduced from dl/dV curves of the smaller 
QDs up to an aerial size A, which corresponds to an average 
radius of r = ^fA/'K. Assuming a circular quantum dot, the 
analytic solution for confined states is given by the zeroes 
of the Bessel function (n = 0,1, 2,...) according to 
Jn{kn • r) = 0. The wave vector kn can be taken from the 
dispersion, which for graphene would read kn = En/idv^) 
with En being the energy of the peak in dl/dV curves. 
Plotting this kn multiplied by r for different QDs shows 
good agreement with the zeroes of the Bessel functions xon 
(xoo = 2.40, xqi = 3.83), if vd = 1.0 • 10^ m/s. This reveals 
that the particular tips used for these QDs preferentially probe 
the graphene states, at least, as it concerns the peak structures 
in dl/dV curves. 

Figure [^b) displays the same data as an E{k) dispersion 
(black squares) using the measured peak energies En and 
kn = ^On/^- The data points fit reasonably with the black 
line of the graphene Dirac cone as measured by ARPES on 
Ir(lll) (311341. However, using the same method except 
that r is determined from the radius of an inner circle which 
completely fits into the island, Altenburg et al find that their 
peaks on average fit better to the Ir(lll) surface state Sq (see 
diamonds in Fig. [^b)), albeit some of the states could also 
correspond to the graphene Dirac cone dispersion (26l. They 
also provided DFT data showing that an s-type tip being 0.48 
A above the graphene is more susceptible to Sq states than to 
states from the Dirac cone. In line, they used an Au tip known 
to be dominated by s-states at E^, while Subramaniam et al 
used a W tip, typically showing stronger contributions from 
d-states. Other groups used different methods to determine 
E{k) for the graphene QDs on Ir(lll). While Phark et 
al (triangles pointing down in Fig. [^b)) used the Fourier 
transformation of the real space dl/dV images resulting in 
an intensity distribution /(^||), which displayed six dominant 
peaks. Angularly averaging /(^| |) results in /(|^| 1 1) and leads 
to a single peak taken as the k value belonging to the energy 
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FIG. 6. (a) Analysis of the QD peak energies in dl/dV curves using a circular model geometry for the graphene QD with radius r determined 
from the measured QD area A by r = such a model exhibits confined states at Jn(A^n • r) = 0 with being the n^^ Bessel function 

and the graphene dispersion En — hv\:,kn assuming ud = 1 .0-10® m/s; the agreement between the experimental symbols and the dashed line 
confirms the Dirac velocity of the confined states to be ud = 1.0-10® m/s O; (b) dispersion of confined states in graphene islands on Ir(l 11): 
Subramaniam et al O and Altenburg et al deduced E{k) from the method described in (a) except that the measured peak energies 
E from dl/dV curves and the kn deduced from kn ■ r — xon with measured QD radius r — and known values of the zeroes of 

the Bessel functions xon are used; Hamalainen et al El compared the dl/dV images with numerically calculated confined wave functions 
using a dispersion of Klein-Gordon type, QD state energies are taken to be the dl/dV images with best agreement in spatial distribution; 
the symbols are then deduced as above using r = with the published QD area A; Phark et al CD used the angularly averaged 

Fourier transformation of the real space data which exhibits a peak in the intensity vs. \k\\ \ representation taken as the k value; the energy 
E corresponds to the voltage V where the dl/dV image is recorded according to E = e ■ V f2 ^ ; the lines are deduced from the ARPES 
data displayed in Fig. [^e) (black line, graphene Dirac cone) and Fig. |^f) (red line, Ir(lll) surface state So with graphene on top); (c) same 
as (b) using the method of Subramaniam et al after intercalating 720 L of O 2 at T ~ 450 K (symbols) in comparison with dispersion from 
graphene Dirac cone according to Fig. [^e) 1^ ; (d) STM image of graphene QD on Ir(lll) after O 2 intercalation; (e)—(g) dl/dV images of 
the graphene QD shown in (d) recorded at the energies E = e -V sls marked. (23 ((a) courtesy of F. Libisch, TU Vienna, (b)-(g) courtesy of 
C. Busse, University of Cologne) 
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FIG. 7. (a) dl/dV image of an interface between Au(l 11) and graphene on Au(l 11) (G/Au(l 11)) as marked with the interface highlighted by 
a dashed line, V = —75 mV, 7=1 nA, Vmod — 3 mV; arrows mark directions of wave guiding visible as a one-dimensional orientation of 
a standing wave fT^I^ : (b) dljdV image on the graphene on Au(lll) with atomic resolution, V = — 20 mV, 7=1 nA, Vmod — 3 mV; 
(c) Fourier transformation (FT) of the dl/dV image in (b) with dashed hexagons marking the Brillouin zone of graphene (smaller hexagon) 
and the larger reciprocal unit cell connecting the spots which originate from the graphene atoms; white arrow is a vector from F to a point 
of the Dirac circle around K, blue arrow marks a vector from F to a state belonging to the Au(lll) surface state; (d)—(f) zoom into the area 
around K at different voltages as marked; (g) ring diameter k deduced from the FT of dljdV images at different V — Ele\ red squares: ring 
around K; dark blue triangles: ring around F; cyan diamonds: ring around F from areas not covered by graphene; lines are linear (red) and 
parabolic (blue,cyan) fits of the equally colored symbols with red dotted lines marking the possible slopes of the red fit line within error bars, 
in (courtesy of M. Fonin, Konstanz University) 


E corresponding to the voltage V = F^/e of the dljdV 
image. The data fit obviously better to the So state, albeit not 
perfectly. Finally, Hamalainen et al. used a comparison of 
apparent lateral shapes in calculation and experiment. They 
visually compared confined squared wave functions in dl/dV 
images and confined squared wave functions from numerical 
calculations using the Klein-Gordon equation. The Dirac 
velocity and, thus, the wave vector k at given E is adapted 
in the calculation until it fits to the experimentally found 
energy, where the calculated state resembles the experimental 
one most strongly (upwards triangles in Fig. [^b)). This 
comparison again results in a better fit of E{k) to So than to 
the Dirac cone of graphene (Note that the claim in the original 
publication was different (Ml). 

It somehow appears that, firstly, the shapes of the confined 
wave functions (used in ED and (Ml) are more strongly 
influenced by So than the peaks in dljdV curves (used in 
a and (Ml)- This points to a longer lifetime of the Dirac 
cone states than of the Ir(III) So state, which makes the 
peaks belonging to the Dirac cone sharper than the peaks 
belonging to So. This suspicion appears to be in line with 
the different energetic widths of So and the Dirac cone found 
in the ARPES data of Fig. [TJe) and Fig. |4];f). However, 
secondly, the relative importance of the two contributions can 
depend on details of the tip which are not under control in the 


STM measurements. 

An easy way to get rid of the So contribution is the 
intercalation of oxygen (750 L O 2 at 450 K) between the 
graphene QD and the Ir(III) (29l . This, firstly, removes So 
as evidenced by ARPES (76l and, secondly, increases the 
graphene-Ir(III) distance. Indeed, Fig. [^c) shows that the 
resulting E{k) dispersion deduced from the peak energies 
and the island sizes, as described for Subramaniam et al. 
above, fits rather nicely to the graphene Dirac cone, which 
is additionally p-doped by the oxygen as also deduced by 
ARPES (Ml . The corresponding confined wave functions, 
shown in Fig. |^e)-(g), exhibit the same regular spatial 
appearance as discussed in Fig. indicating that the soft 
confinement driven by the edge-Ir(III) interaction is still 
present even after oxygen intercalation. 

While this line of probing of graphene QD states has not 
been pursued further so far, an interesting way to distinguish 
the contributions from Dirac cone states and interfering 
surface states has been published (Z7l . In order to decouple 
the graphene QDs from the Ir(III), the authors intercalated 
5-10 nm of Au(lll) between the QDs and Ir(III) by Au 
evaporation at 300 K. Indeed, the graphene islands on top 
of Au(lII) could now easily be moved by the forces of 
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the tip of the STM indicating that strong chemical bonds 
as between Ir(lll) and the graphene edges are not present 
on Au(lll). Nevertheless, the dominating Au(lll) surface 
state, being located around T and intersecting Ey CTl |78l, 
still interferes with the Dirac cone states. Figure |7Ja) shows 
a dljdV image with an interface (dashed line) between 
uncovered Au(lll) areas (left) and Au(lll) areas covered 
with graphene (right). Both areas show very similar standing 
wave patterns. Interestingly, the wave guiding along the 
herringbone reconstruction (arrows) 1791 . recently explained 
as a diffraction effect l80l . is visible, too, on both sides of the 
interface. Thus, obviously the Au(lll) surface state, known 
to have a parabolic dispersion as E{k) = j 
with an origin at E^ — Ey = —480 meV and an effective mass 
of = 0.25 • me (jn^. electron mass) 17711781 . is visible by 
tunneling on the graphene revealing that it is penetrating into 
the graphene like the Sq state of Ir(l 11 ). 

However, by recording dl/dV images with atomic resolution 
as shown in Fig. |^b), it gets clear that the graphene Dirac 
cone contributes to the observed standing wave patterns 
as well. The Fourier transformation (FT) of the real-space 
dl/dV image (Fig. [^c)) shows a central circle compatible 
with the Au(lll) surface state and, in addition, circles which 
correspond to intervalley scattering between K and K’, i.e. 
between the Dirac cones of graphene. These circles show 
an energy dependence of its diameter as expected for the 
graphene Dirac cone (Fig. [7Jd)-(f)), i.e., while the inner 
circle belonging to the Au(lll) surface state gets larger with 
increasing energy, the outer circle gets smaller. Plotting the 
diameter of the circles as a function of energy (Fig. [TJg)) 
allows to deduce the dispersion of the contributing bands. 
The authors find a parabolic dispersion for the inner circle 
with an effective mass fitted to be m* = 0.26 ± 0.02 • me 
as found identically on the uncovered Au(lll) areas. This 
identifies this inner ring as being caused by the Au(lll) 
surface state penetrating into the graphene. The fitted band 
onsets, Eq — Ey = —390 meV on clean Au(lll) and 
Eq — Ey = —300 meV on graphene covered Au(lll), are 
higher than on the Au(lll) surface of a single crystal, which 
is probably due to confinement effects of the thin Au(lll) 
layer on the Ir(lll) and charge transfer from the graphene, 
respectively. The diameter of the outer rings reveals a 
linear dispersion E = hv^k + E^ with fit parameters 
vd = (1.1 ± 0.2) . 10® m/s and Eu = 0.24 ± 0.05 eV 
identifying it as belonging to the graphene Dirac cone. 
Thus, in principle, the intensity ratio between the inner and 
the outer rings provides access to the relative contribution 
of Dirac cone states and surface states of the support to 
the dl/dV images of the graphene QDs. It also allows to 
discriminate their spatial shapes by an adequately filtered 
back transformation of the FT data, an interesting experiment 
still to be done. 


C. Edge states 

Coming back to the graphene QDs on Ir(lll), we discuss 
its edge properties. As shown in Fig. [Jb), the edge is com¬ 
pletely zigzag type. This might be regarded as favorable con¬ 
sidering that the armchair edge is much more stable for freely 
suspended graphene |[40| [811 and that graphene zigzag edges 
are predicted to exhibit an edge state becoming ferromagnetic 
if Ey — Ey (3911121 • However, this edge state is suppressed 
on Ir(l 11 ) as will be explained in the following. 

DFT calculations of a graphene ribbon on top of Ir(l 11) show 
that the strong interaction between the Ir d;^ 2 -orbitals and the 
Pz orbitals of graphene lead to a preferential on-top position 
of the C-atoms at the edge with respect to Ir atoms (see re¬ 
laxed geometry in Fig. [^a), left area). The bond induced 
downwards bending of the C atoms by about 1.6 A leads to a 
sp^ like bonding configuration of the C atom with two bonds 
to neighboring C’s and one bond to the neighboring Ir d ^2 
orbital. The remaining sp^ bond is most likely saturated by 
hydrogen which originates from the preparation via C 2 H 4 . 
The latter has been verified by comparison between mea¬ 
sured STM images at different bias and simulated STM im¬ 
ages based on DFT (3^ . Moreover, the single-H terminated 
C-atoms at the edge are energetically favorable with respect 
to double-hydrogenated C-atoms and unsaturated C atoms on 
top of Ir(lll) (361 . Importantly, the preferential on-top po¬ 
sition of the C-atoms causes strain within the graphene edge 
region due to an effective lattice mismatch between graphene 
and Ir(lll) of about 10 %. This, in turn, causes strain relax¬ 
ation areas, where the graphene edge is more detached from 
the Ir(l 11 ) being close to its natural distance in the interior of 
extended graphene of 3.4 zb 0.2 A (831 (see right area of Fig. 
[^a)). These more detached areas mimic a sp^ bonding con¬ 
figuration with two bonds connected to neighboring C atoms 
and the remaining one being saturated by H, while the edge 
Pz orbital remains unsaturated. 

Surprisingly, the DFT calculations did not find indications of 
the predicted edge state being a hallmark of the zigzag edge. 
This is shown by the projected density of states (PDOS) of 
different atoms around the edge (Fig. [^b) and (c)), which 
do not show very pronounced features in the area around Ey, 
i.e. Sit E = 0 eV. While the C atoms in the sp^ areas show a 
suppressed PDOS with respect to the surroundings, the atoms 
from the sp^ area show a triple peak at the very last atom, 
which is, however, not larger than other peaks of the PDOS. 
For comparison, the red curve in Fig. [^d) shows a typical 
edge state which is a factor of 5 — 10 higher than all surround¬ 
ing peaks. In line, the dl/dV curves recorded at the edge of 
the graphene QDs on Ir(l 11) do not show any spectral feature 
around Ey as marked in Fig. [^e). This result is found for 
several tens of graphene QDs using a similar number of dif¬ 
ferent microtips. 

The absence of the edge state is surprising, since it is generally 
believed that the edge is a relatively robust feature of graphene 
zigzag edges. It originates from the fact that the introduction 
of radical tt- bonds on every third edge atom allows the forma¬ 
tion of Clar’s sextets (six double bonds in one benzoid ring) 
through the whole interior of the graphene bulk, such that sub- 
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FIG. 8. (a) Atomic arrangement of a graphene nanoribbon on the Ir(l 11) substrate (blue balls: Ir atoms, black balls: C atoms of graphene bulk, 
red balls: graphene edge atoms, pink balls: H atoms); the arrangement is deduced from the relaxation within a DFT calculation; for the sake 
of clarity, only the top Ir layer is shown, while four layers of Ir are used in the calculations; different atoms are numbered and the bonding 
configuration (sp^ or sp^) of several C atoms at the edge is marked; (b)-(c) projected density of states (PDOS) of particular atomic spheres; 
(b) atoms behind the sp^ bond edge atom as indicated by the red line; rows are numbered away from the edge; (c) same as (b) for the atoms 
behind the sp^-type edge atom marked by a red line, too; (d) PDOS of C edge atoms (as marked by lines) of a freely suspended graphene 
nanoribbon as sketched in the inset with partial twofold H-termination (colored areas) and partial one-fold H-termination (cyan H atoms at 
black C atoms); the red and blue colors in the inset mark the oppositely oriented spin polarizations of the atoms with the size of the colored 
ball marking the strength of spin polarization; (e) dl/dV curves recorded at the positions marked in the STM images of a graphene QD in the 
inset, /stab = 0.3 nA, I/gtab = 0.3 V, Vmod = 5 mV; the area around Er>, where an edge state would be expected, is marked. ((a)—(d) 
courtesy of R. Mazzarello, RWTH Aachen University) 


stantial kinetic energy by the corresponding delocalization of 
electrons within the pz orbitals of the ring can be gained 1401 . 
The radicals at the edge then form the flat band at be¬ 
ing prone to a ferromagnetic instability, if E^ = E^ 1^ . 
Thus, at first hand, the soft confinement found in the previous 
paragraph for graphene QDs on Ir(III) would just move the 
necessity of these radicals to the interior of the graphene and, 
thus, would lead to an edge state away from the physical edge 
of the graphene QD. Such an edge state, however, is also not 
observed neither in DFT nor in STS. 

There are two possibilities to explain the disappearance of 
the edge state. Firstly, the bonding along the edge is not 
translationally invariant on the atomic scale, but resembles a 
zz(2222222III) configuration in the nomenclature by Wass- 
mann et al. 1401 where zz symbolizes the zigzag edge and 
the numbers describe the number of saturated bonds in the 
periodically repeated sequence of ten edge atoms. Our DFT 
calculations of a freely suspended graphene ribbon with such 
an edge termination, represented by H atoms, reveal that this 
termination, indeed, partly destroys the edge state, i.e. the 
radicals are only required in the area of the double bonds in 


order to allow the formation of Clar’s sextets 1361 . This can 
be rationalized by the fact that a zigzag edge length of three 
unit cells can host Clar’s sextets in the whole interior of the 
graphene ll40l . However, for an extended edge, this requires 
that the surrounding edge areas host ionized pz bonds to adapt 
Clar’s sextets in the rest of the graphene bulk. Secondly, the 
presence of the substrate can prefer certain locations of dou¬ 
ble bonds by the interaction with the dz 2 orbital of the Ir(111), 
such that the gain in kinetic energy for Clar’s sextets is over¬ 
compensated by the exchange energy between Ir dz^ orbitals 
and the C pz orbitals. Such an effect should be directly visible 
in STM images which would show these preferred bonds as 
occupied states close to E^, while the unfavorable bonds are 
visible as unoccupied states mni. 

Indeed, differences between brightly appearing bonds in 
occupied and empty state STM images (negative and positive 
bias, respectively) are visible in Fig. While the bonds per¬ 
pendicular to the zigzag edge are preferentially empty, other 
bonds more parallel to the zigzag edge appear to be occu¬ 
pied. Thus, we conclude that the lattice mismatch between 
Ir(III) and graphene, leading to the zz(2222222III) config- 
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FIG. 9. (a)—(d) STM-image of edge of graphene QD on Ir(l 11) with 
atomic resolution; (a), (b) 1/ = —0.5 V, / = 5 nA; (c), (d) V = 
—0.5 V, / = 5 nA; white lines mark preferential bond orientations 
within the edge regions of the images. 


uration and to the inequivalence between C-C bonds on the 
substrate in combination with the ( 1^2 orbital at Ep, leading 
to a relatively strong interaction with the graphene orbital, 
suppresses the edge state on Ir(lll) 1361 . This is favorable 
for graphene QDs to be used as spin qubits nni, since the 
edge magnetism could interact with the qubit spin, thus, of¬ 
fering a prominent channel for decoherence. However, it is 
not straightforward to transfer this concept of edge state sup¬ 
pression to graphene QDs on insulating substrates. 


TABLE I. DFT results regarding the presence of magnetic edge states 
at graphene zigzag edges on different metallic substrates with differ¬ 
ent terminations Ea. 


substrate/ edge state magnetic 
termination 

reason 

Ir(lll)/H 

No 

- 

hybrid with Ir 5d_22 

Ir(lll)/bare 

No 

- 

hybrid with Ir 5d^2 

Au(lll)/H 

Yes 

Yes 

- 

Au(lll)/bare 

Yes 

No 

charge transfer to Au 

Ag(lll)/H 

Yes 

No 

charge transfer to Ag 

Ag(lll)/bare 

Yes 

No 

charge transfer to Ag 

Cu(lll)/H 

Yes 

No 

charge transfer to Cu 

Cu(lll)/bare 

Yes 

No 

charge transfer to Cu 


In order to demonstrate that the suppression of the edge 
state is rather the exception than the rule, we performed cal¬ 
culations of graphene nanoribbons on different close-packed 
metallic surfaces such as Ag(lll), Cu(lll) and Au(lll), ei¬ 
ther with single H-termination of the edge or without termina¬ 
tion (53). The results in combination with the ones on Ir(l 11) 
are shown in table All zigzag edges except the ones on 
Ir(lll) exhibit an edge state, most likely, since none of these 


surfaces has a surface state close to Ejj. Thus, the pres¬ 
ence of this particular surface state pointing directly to the 
graphene orbital and being partly occupied, since at Ep, 
is decisive for the suppression of the edge state as analyzed 
above. 

Albeit most of these substrates leave the edge state intact, the 
charge transfer from the edge states to the substrate mostly 
results in an empty edge state, which is not prone to a fer¬ 
romagnetic instability, as consistently found in the DFT cal¬ 
culations (531. Only a Au(lll) substrate, which exhibits the 
weakest interaction and only after an additional H-termination 
of the graphene edges, which further reduces the interaction 
with the Au(l 11) (1^2 surface states being located 2 eV below 
Ep, leaves the charge transfer small enough, such that a ferro¬ 
magnetic instability occurs according to DFT (361. In line, the 
published indirect evidence for a ferromagnetic edge state (44l 
was observed on Au(l 11), however, without a controlled edge 
termination. Subsequent hydrogenation of the corresponding 
graphene nanoribbon edges has also been achieved (55), how¬ 
ever, so far, without an investigation of the magnetic proper¬ 
ties of the edge. 

Also in line, the edge state of the graphene QDs on 


-60 -40 -20 0 20 40 60 

Bias (mV) 

FIG. 10. (a) STM-image of the lateral zigzag interface between 
graphene and BN as marked on Au(lll)/Ir(lll) showing the edge 
state as a bright area, V = —1 V, I = 3 nA; (b) dl/dV curves 
recorded on graphene, BN, and the lateral interface as marked; (c) 
stick-and-ball model of a freely suspended lateral graphene-BN in¬ 
terface with C-atoms (grey), B-atoms (orange), and N-atoms (blue); 
(d) DFT based simulation of the spatial distribution of the tunnel¬ 
ing current of the assembly of (c) at 1/ = —0.2 V with positions of 
C-atoms (red), B-atoms (orange), and N-atoms (blue) overlaid. (601 
(courtesy of P. Liljeroth, Aalto University) 

Ir(lll) can be recovered by intercalating Au(lll). There¬ 
fore, Drost et al. io), firstly, terminated the graphene zigzag 
edges by BN, which can be grown on Ir(lll) using borazine 
((HBNH) 3 ) (85l and which seamlessly contacts the previ¬ 
ously prepared graphene QD edges along the zigzag direc¬ 
tions. Secondly, they intercalated Au such that the interaction 
with the substrate is additionally reduced by the removal of 
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the dz 2 surface state of Ir(lll). Figure [T^a) shows a STM 
image of the lateral interface between graphene and BN on the 
Au(lll) showing a state at the interface which slightly varies 
in perpendicular size probably due to a remaining influence of 
the substrate. The corresponding state could be identified by 
dljdV curves (FIG.p^b)). It is found to be within ±5 meV 
around in accordance with DFT calculations tbOl . It is not 
clear so far, if the state is ferromagnetic, which can be probed, 
e.g., by spin-polarized STM 1861 . However, the lateral distri¬ 
bution of the state can be nicely reproduced by DFT ignoring 
the substrate as shown in Fig. [^c) and (d), if the energet¬ 
ically preferential B-termination of the BN-graphene zigzag 
interface is chosen 1601. 

Thus, an edge state can indeed be recovered from the favor¬ 
able production of zigzag edges on Ir(l 11) being located close 
to Ey. It might exhibit the long-sought edge magnetism di¬ 
rectly. This would be not favorable for the spin qubit appli¬ 
cations, but other spintronic applications as, e.g., spin Alters 
relying on this edge state have been proposed l56l . 


II. ANODICALLY OXIDIZED QUANTUM DOTS 

The described STM experiments on graphene QDs de¬ 
posited on metals provide an unprecedented spatial resolu¬ 
tion concerning the electronic properties of the QDs, but they 
are difficult to be directly correlated to transport properties 
of graphene QDs, which are typically measured on insulat¬ 
ing substrates and which are defined in a top-down approach 
by etching and subsequent contacting of graphene flake struc¬ 
tures using e-beam lithography ca. It is believed that the 
difficulties to control the QD transport properties in a similar 
way as in GaAs Ell is mostly due to the uncontrolled edge 
geometry and chemistry |[2T] |88l . This is distinct to the QDs 
studied in UHV by STM, where these two parameters are well 
controlled. 

In order to locally probe transport QDs, scanning gate mi¬ 
croscopy has been used, which, however, suffers from a rather 
limited spatial resolution of about 100 nm |[6T][62l. STS has a 
much higher spatial resolution, but requires that the QD is on 
the surface, since the tunneling current decays by an order of 
magnitude, if the tip-QD distance is increased by merely 1 A 
1891 . Moreover, the surface of the QD has to be free of insu¬ 
lating resist, which otherwise would be charged by the electric 
held of the tip. In turn, this charging of resist can change the 
local potential below the tip and, thus, the LDOS contributing 
to the tunneling current l90l . Consequently, it is mandatory 
to structure and contact the graphene nanostructures without 
resist. 

For contacting, one can either use shadow mask evaporation 
or contacting by microsoldering 1^ . We used the second 
method l68l . which basically melts a fine wire of In onto the 
graphene sample held at about 155^ C. The In wire is pulled 
out of an In droplet by a typical STM tip made out of W. The 
positioning of this In wire onto the sample is afterwards done 
with micromanipulators under an optical microscope. This 
leads to a multitude of In contacts placed with a precision of 
1-10 jim as shown in Fig. [^a). In order to structure the 


graphene into a QD, we used local anodic oxidation (LAO) 
by the tip of an AFM 16311 . This technique has successfully 
been used previously to produce nanostructures in GaAs 1921 . 
It also works on graphene l64l . where, e.g., quantum dots l65l 
or constrictions l93l made by LAO have been probed by trans¬ 
port spectroscopy. 

The basic technique is to move the cantilever of an AFM 
across the surface with a velocity vnp ^100 nm/s while ap¬ 
plying a negative tip bias Vap. A controlled humidity leads to 
an adhesive droplet of water between tip and sample, where 
the electric held between tip and sample dissociates the H 2 O 
and guides the oxygen radical to the sample surface, such that 
an oxidation of the graphene can take place. The C-oxidation 
on graphene can be either complete leading to cuts in the 
graphene l64l or incomplete, most likely leading to partially 
oxidized graphene nanoislands, which provide insulating lay¬ 
ers, too (Sl|67l. These partially oxidized areas appear el¬ 
evated with respect to the intact graphene in tapping mode 
AFM images. They are electrically insulating up to voltages 
of about 2 V 1661 . Typically, cuts are surrounded by incom¬ 
pletely oxidized areas. Reducing the tip voltage, in order to 
selectively produce the oxidized nanoislands, leads to the lat¬ 
erally smallest insulating barriers. 

As shown in Fig. [mb), we have produced a graphene QD 
with a diameter of about 60 nm using LAO. This QD contains 
source (S) and drain (D) contacts and a number of plunger 
gates (pg) in order to tune the electrostatic potential of the QD 
and its tunneling barriers. The QD is placed on a Si/Si02 sub¬ 
strate such that the highly doped Si can act as a backgate. We 
could show that the etching process of graphene by LAO also 
works on BN supports. 

Figure [TTJc) shows the measured 2-point conductance G = 
Isb/Vsb (Vsd' voltage applied between S and D, Jsd- cur¬ 
rent measured in the circuit including S, D, and QD) through 
the QD as a function of the backgate voltage I4ackgate re¬ 
vealing a Vbackgate area, where the conductance in between 
peaks drops towards zero as typical for Coulomb blockade 
ISTll . Zooming into this fbackgate area exhibits that the dif¬ 
ferential conductance dIsB /<^Vsd indeed drops significantly 
between peaks of full width at half maximum (FWHM) of 
AVbackgate = 15 mV. Plotting the differential conductance as 
a function of fbackgate and Vsd reveals the typical Coulomb 
blockade structures (Fig. [^e)) known for QDs |[T4l[T5ll87l . 
However, the patterns in Fig. [^e) appear like overlapped 
Coulomb structures from QDs connected in series. Only a 
few of the diamonds like the one marked by a white dot have 
the characteristic shape of an isolated charging event and can 
be used to determine, e.g., the lever arm of the back-gate volt¬ 
age to the QD to be a = 0.067. Most likely, our structure still 
suffers from additional confinement areas within the source 
and drain regions, which produce additional quantum dots in 
series, as partly observed in etched graphene QDs, too 1941 . 
The size of the marked Coulomb diamond fits reasonably to 
the QD area of Fig. pT]^b). Taking a QD area of A = 60 x 200 
nm^, a thickness of the Si02 of d = 300 nm and a dielectric 
constant of e = 4, we get an estimated back-gate-QD capaci- 
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FIG. 11. (a) Optical micrograph of a graphene flake on Si/Si02 contacted by seven In lines microsoldered onto the graphene flake ll^ : (b) 
AFM tapping-mode image (phase image) of a graphene QD (QD) on Si/Si02 with source (S) and drain (D) contacts as well as a number of 
plunger gates; the plunger gate used in (f) is marked as pg; the oxidized lines (bright lines) are produced by LAO in AFM contact mode at 
a controlled humidity of 60 %, Vtip = —8.5 V, utip = 200 nm/sec; (c) conductance as a function of Vbackgate at Vsd = 1.5 mV and 
Vpg = — 20 mV; all other contacts are grounded, T = 0.3 K; transport gap where G partly drops to zero is marked; (d) d/so/o^VsB as a 
function of Vbackgate within the transport gap region; an oscillating Vsd with amplitude VsD^ax = 0-1 is used for lock-in detection; 
all other contacts are grounded, T = 1.2 K; (e) logarithmic representation of ti/sD/dVsD as a function of Vsd and Vbackgate with all other 
contacts grounded; an oscillating voltage with VsD^ax = 0-1 is added to Vsd for lock-in detection, T = 1.3 K; the white dot marks 
a well-developed Coulomb diamond discussed in the text; (f) logarithmic representation of d/sD /dVsn as a function of Vsd and Vpg with 
Vbackgate = 14.8 V and all other contacts grounded; an oscillating voltage with VsD^ax = 0-03 mV is added to Vsd for lock-in detection, 
T = 0.3K. 


tance Cbg using the plate capacitor model: 

Cbg = = 1-4 aF. (2) 

Dividing this by the lever arm a results in a total capacitance 
of Ctot = a = 21 aF such that we get a charging voltage 
of AU = e I Clot = 7.6 mV in reasonable correspondence 
to the 4 mV found for the marked Coulomb diamond in Fig. 
pT]^e). The remaining discrepancy is probably due to the in¬ 
adequate plate capacitor model, since the size of the QD is 
smaller than its distance to the backgate leading, e.g., to fringe 
field amplification ||95]| . 

Fig. [TTJf) shows that the graphene QD can also be tuned 
by the plunger gate with a much better lever arm close to 
a = 0.2, which is a typical value of a as also found for 
etched graphene QDs ifld^ITTl 1941 . Using the charging volt¬ 
age of the plunger gate to estimate its capacitance, we get 
Cpg = e/AUpg 3.5 aF (AU^g 45 mV) such that the 
six gates as well as source and drain can easily account for the 
lever arm of the back gate. 

Thus, the anodic oxidation can produce graphene QDs tun¬ 
able by backgate and sidegates, where the latter exhibit and 
extremely strong tuning ability. Moreover, the QDs are pro¬ 


duced without any resist, such that they can be probed by 
transport spectroscopy and scanning tunneling spectroscopy 
at the same time. A scanning tunneling microscope having 
five transport contacts for the sample, operating in UHV at 
0.4 K, and in magnetic fields up to 14 T is available (961, such 
that the mentioned experiments to map the quantum dot states 
responsible for transport are within reach. However, it would 
be favorable to probe quantum dots with isolated Coulomb 
diamonds only, using, e.g., a BN support (23. 

III. CONCLUSION 

We have shown, that scanning tunneling spectroscopy of 
graphene nano-islands on metallic substrates, in particular on 
Ir(l 11) and Au(l 11), can probe the squared wave functions of 
graphene QDs, which are, however, intermixed with metallic 
surface states penetrating into the graphene layer. Effective 
ways to decouple the surface states are either, physically, by 
intercalating oxygen between graphene and the metal or, an¬ 
alytically, by performing a Fourier transformation provided 
that the real space images exhibit atomic resolution, such that 
the contributions from graphene states located at K and K’ 





























13 


and the metallic surface states located at F get distinct in 
Fourier space. First investigations of the graphene QD states 
on Ir(lll) showed that the soft confinement induced by the 
interaction between the graphene p^ orbital and the Ir d ^2 or¬ 
bital at leads to rather regular confined states, which are 
not prone to a detrimental K-K’ mixing by the edge of the QD. 
Moreover, this interaction in combination with the 10% lattice 
mismatch between Ir(lll) and graphene suppresses the edge 
state, otherwise, typical for zigzag edges of graphene, such 
that it can not interfere with the spin properties of the confined 
state. In contrast, the edge state is predicted to be present, e.g., 
on Cu(lll), Ag(lll) and Au(lll) supports and it can be re¬ 
covered by laterally interfacing the zigzag edges on Ir(lll) 
with BN and subsequently intercalating Au. Thus, a tunabil- 
ity of the graphene QDs in a favorable or unfavorable way by 
the interaction with the substrate and/or the lateral interfaces 
has been demonstrated, which might be a helpful guide to tune 
quantum dot properties also for applications. 

In order to go beyond these model studies of graphene QDs 
under well controlled conditions, we finally presented an ap¬ 
proach to combine transport studies and scanning tunneling 
spectroscopy directly. Therefore, we use local anodic oxi¬ 
dation for nanostructuring and subsequent micro soldering for 
contacting, such that the nanostructures are free of resist as 


required for the scanning probe studies. It will be interest¬ 
ing to see in how far this novel approach will ease an opti¬ 
mization of graphene QDs, e.g., with respect to applications 
as qubits, where graphene has large fundamental advantages 
with respect to other materials as GaAs, but experimentally 
still lags significantly behind. 
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